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Mixing property and pseudo random 
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Nihon University 

Abstract: We will give a summary about the relations between the spectra of 
the Perron-Frobenius operator and pseudo random sequences for 1-dimensional 
cases. 



There are many difficulties to construct general theory of higher-dimensional cases. 
We will give several examples for these cases. 

1. Perron Frobenius operator on one-dimensional dynamical systems 

Let / = [0, 1 } and F be a piecewise monotonic and C 2 transformations on it. We 
assume F is expanding in the following sense: 

£ = liminfessinf -log\F n '(x)\ > 0. 

n—>oo x£i n 

Definition 1.1. An operator P: L 1 — > L 1 defined by 

V- F{y)=x 

is called the Perron-Frobenius operator associated with F. 

It is well known that the spectra of the Perron-Frobenius operator P associated 
with the transformation F determine the ergodic properties of the dynamical system 
(I,F). Roughly speaking: 

1 . The dimension of the eigenspace associated with eigenvalue 1 equals the num- 
ber of the ergodic components. 

2. The cigcnfunction p of P associated with the eigenvalue 1 for which p > and 
/ pdx = 1 becomes the density function of the invariant probability measure. 

3. If 1 is simple and there exists no eigenvalue modulus 1 except 1, then the 
dynamical system is mixing. 

4. The second greatest eigenvalue determines the decay rate of correlation. 

To be more precise, the Perron-Frobenius operator satisfies: 

1. P is contracting, and 1 is an eigenvalue. Hence, its spectrum radius equals 1. 

2. P is positive. Hence there exists a nonnegative eigenfucntion associated with 
the eigenvalue 1. 
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3. On th unit circle, the eigenvalues of P coincide with the eigenvalues of the 
unitary operator Uf(x) = f(F(x)). 

4. The essential spectrum radius of P: L 1 — > L 1 equals 1, that is, it is of no use 
to consider decay rate for the functions in L 1 . 

5. When we restrict P to BV the set of functions with bounded variation, then 
its essential spectrum radius equals e - ^. Hence, the decay rate for a function 
with bounded variation is at most e - ^. 

Hereafter, we only consider P restricted to BV. Since P is not compact, it is not 
easy to calculate its spectra. Hofbauer and Keller (Q) proved that they coincide 
with the singularities of the zeta function 



CO) = exp 



oo n 
n=l p: F n (p)=p 



The zeta function has the radius of convergence 1, and it has a meromorphic 
extension to But it is not easy to study the singularities in \z\ > 1. 

Since F is piecewise monotone, there exists a finite set A which we call alphabet 
with the following properties. 

1. For each a £ A, an interval (a) corresponds, and {(a)} ae _4 forms a partition 
of/, 

2. F is monotone on each (a), and it has a C 2 extension to the closure of (a). 



We define 

+ 1 if F is monotone increasing in (a) 
— 1 otherwise. 



sgna = 



We call F Markov if there exists an alphabet A, such that, if F({a))° (~l {b) ^ then 
F((a)) D (b)°, where J° denotes the inner of a set J. 

We will summarize notations which we use in this article. 

1. A sequence of symbols a± • ■ • a n (<H £ A) is called a word and 

(a) |io| = n (the length of a word), 

(b) (w) — n"=i F~ l+1 ((a,)) (the interval corresponding to a word), 

(c) We denote the empty word by e, and define |e| = 0, sgne = +1 and 
(e} = /, 

(d) we call w is admissible if (w) ^ 0. We denote the set of admissible words 
with length n by W„ and W = U£L W„, 

(e) sgnw = nr=i s S na i- 

2. For a point x £ /, af af • • ■ is called the expansion of x which is defined by 

F n (x) £ (a n+1 ). 

We usually identify x and its expansion. 

3. For a sequence of symbols s = a\a-2 ■ ■ ■ (a, £ A), 

(a) s[n,m] — a n a n+ i ■ ■ ■ a m (n < m), 

(b) s[n] = a n . 
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4. For a word w = a\ ■ ■ ■ o„ and x £ /, is an infinite sequence of symbols 
defined by 

wx = a\ ■ ■ ■ a n a\a^ ■ ■ ■ . 

We call wx exists if there exists a point y whose expansion equals wx. 

We have a natural order on A, that is, a < b when x < y for any x € (a) and 
y 6 {b}. Then we define an order on words by 

1. w < w' if |w| < \w'\, 

2. w — a\ ■ ■ ■ a n , w' = b\ ■ ■ ■ b n and if = 6,; for k < i < n and au < bk, then 

w < w' if sgnafc+i • • • a„ = +1, 
w > w' otherwise. 



2. Pseudo random sequences 

We call a sequence {x n }^° =1 uniformly distributed if for any interval J C / 

#{x n eJ:n<N} 

J™ M = J ' 

where \J\ is the Lebesgue measure of an interval J. 

Uniformly distributed sequece {x n }^ =1 is called of low discrepancy if 



Dn = sup 

J 



#{x n eJ:n <N} _ 



N 



O 



logiV 
N 



It is well known that this is best possible and low discrepancy sequences play 
important role in numerical integrations. 

We can construct a low discrepancy sequence using dynamical system. 

Definition 2.1. The set {wx} we w for which wx exists and arranged in the order 
of words is called a van der Corput sequence defined by F . 

When F(x) — 2x (mod 1), A has two elements so we express A = {0, 1}, and 
every words are admissible. Thus our van der Corput sequence is 

x, Ox, lx, OOx, lOx, Olx, llx, 000a;, . . . , 

and original van der Corput sequence is the case when x = i . 
Our first goal is the following theorem: 

Theorem 2.2 ([4]). Let F be a Markov and topologically transitive transformation 
with the same slope \F'\ = (3 > 1. Then, for any x £ I, the van der Corput 
sequence is of low discrepancy if and only if there exists no eigenvalue of the Perron- 
Frobenius operator in the annuls < \z\ < 1 except z = 1. 

For Markov /^-transformations, Ninomiya (|2|) proved the necessary and suffi- 
cient condition that the van der Corput sequence for x = is of low discrepancy by 
direct calculation. He extended the results to non Markov /3-transformations ( 10] ). 

Rough sketch of the proof is the following. For any interval J 

p"ij(x) = Yl ^Ay^'ivT 1 

= (3- n U(wx). 

\w\—n 
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Note that 53|«)|=n 1j( wx ) is the number of hits to J in the van der Corput sequence 
corresponding to words with length n. On the other hand, as a rough expression, 

P n lj(x) = \J\p(x) + r) n , 

where 77 is the second greatest eigenvalue in modulus. Thus the number of wx which 
hits to J with \w\ = n equals /3 n \J\p(x) + (/3ri) n . The discrepancy depends on the 
term (/3ry) n , and we already knew that \rj\ > Thus r\ = /3 _1 is the best possible 
case, and at that time, the number of wx which hit to J with \w\ < n equals 

(3 n+1 - 1 
P_ x \J\p(x)+n. 

The number N of admissible wx with | w| < n is of order /?". This says the discrep- 
ancy is of order lo ff N . 



3. Predholm matrix 

To make the proof of Theorem 12 . 21 rigorous, we define generating functions. Let F 
be a piecewise linear Markov transformation with the same slope [3 > 1. Let for an 
interval J 

00 

s J {z,x) = z n P n lj{x). 

71 = 

Let s(z, x) be a vector with coefficients s^{z^x) (a G .4). Then we get a renewal 
equation of the form: 

S (z 1 x) = (I-^(z))- 1 x(x), 

where <&(z) is a A x A matrix and xi x ) is a dimensional vector whose coefficients 
equal: 



*(2)a,. 




if J F«a))D<6)°, 
otherwise, 

X(x)a = 1(a) (») < 

We call the matrix $(z) the Fredholm matrix associated with F. 
Theorem 3.1. 

det(J - $(«)) ' 



This theorem says that the eigenvalues of P is determined by the zeros of det(7 — 
$(z)). We can express s J (z,x) using the Fredholm matrix and a row vector whose 
coefficients are the functions with the radius of convergence j3 (0|)- 

From these, we can get the rigorous proof of Theorem 12. 21 

We can extend the above discussion to more general cases for 1-dimensional 
transformations. 

We use the signed symbolic dynamics defined by the orbits of the endpoints of 
(a) (aE A). 
Let for x 6 I 

x + = lim a\a^ • • • , 
x~ = lim a\a\ ■ ■ ■ . 

ylx 
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For an interval J, we denote (sup J) + and (inf J)~ by J + and J~, respectively. 
Especially, for a £ A, we denote by a + and a~ instead of (a) + and (a) - . 

Let .A be the set of a a (a £ A,cr = ±). 

We also define 

oc 

S y"(z,x) = zU P~ n E cj{y%wx)8[{w[l}) D (c%),36wx], 
n=o i«ew„ 

where 9 is the shift and 



5\L] 



1 if L is true, 
if L is false, 



1-2 li 2/<<7 a; I 



x < c y 



x < y a = +, 
x > y a = — . 



Then 

s J (z,x) = s J (z,x) + s J (z,x). 
We can also get the similar equation 



s(z,x) = (/-*(«)) 1 x{z,x), 



where 



s(z 1 x) = (s a (z,x))aeA,tr=±, 
X( z , x ) = (X^ (z,x))aeA,a=±, 

'a{y a ,x) if 0y a = (a 2 ') <TS s na i 

k XT=o z n (3- n <j{6 n y a , x) otherwise, 

®{Z) = (0a-,6-(z))a,6e^, ( j,r = ±, 

!+^-V2 if 6 r < CT 6^, and 6y a = (af) <TS s nQ i 
-z[3~ 1 /2 ifb T > a 9y a , and fly* = (af) <TS s nQ i 
£^ =1 z"/rV(6>"?/ cr ,6 T )(sgn 2 / T [l,n- 1]) otherwise. 

Moreover, Theorem 13.11 also holds (cf. 0). This leads to the general theory con- 
cerning the discrepancy of sequences generated by one-dimensional piecewise linear 
transformations. 

Definition 3.2. We call an endpoint of (a) a Markov endpoint if the image of this 
point by F n for some n > lcoincides with some endpoint of (b) (b £ A). 

Theorem 3.3 ([2]). Let us denote by k the number of non-Markov endpoints of 
F. 

1. Let CnP~ n be the n-th coefficient of (1 — z)£(z). Assume that is bounded, 
then 

'(logA0 fe+2> 



D N =0 



N 
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2. Let <I>ii (z) be the minor of the Fredholm matrix <&(z) corresponding to Markov 
endpoints. Assume that det(7 — $n(z)) ^ in \z\ < f3 and there exists no 
sigularity of C( z ) i> n \ z \ < P except 1. Then the discrepancy satisfies 



4. Higher-dimensional cases 

We want to extend the result to higher-dimensional cases. 

Let I = [0, l] d . We will call a product of intervals also by intervals. We can 
define the Perron-Frobenius operator P for higher-dimensional cases. Here we use 
Jacobian detF instead of derivative F'. For Markov cases, we can find the spectra 
of the Perron-Frobenius operator by the Fredholm matrix by almost the same 
discussion. Even for non Markov piecewise linear cases, instead of signed symbolic 
dynamics using the idea of screens we can construct Fredholm matrix and we can 
determine the spectra of the Perron-Frobenius operator (@). 

We consider a space of functions for which there exists a decomposition / = 
J2 W @wl(w) such that for any < r < 1 



This is a slight extention of the space of functions with bounded variation in 1- 
dimensional cases. Here, note that, though the corresponding zeta function has 
meromorphic extention to \z\ < e*, even for the simplest Bernoulli transformation 
F(x,y) = (2x,2y) (mod 1), the essential spectrum radius is greater than e - '. 
Here is a question. 

Are there any transformation whose spectrum radius equals e~^? 

If such a transformation for which the Jacobian is constant exists, we will be able 
to construct van der Corput sequences of low discrepancy. Here we call a sequence 
of low discrepancy if 



where sup is taken over all intervals J C /. 

It is proved that this is best possible for d = 2, and this will be also true for 
d > 3. 

We have not yet the answer, but we have constructed low discrepancy sequences 
for d — 2 using dynamical system. Also for d = 3 we have constructed low discrep- 
ancy sequences using a sequence of transformations F%, . . . but we have not yet 
constructed a transformation for which F n — F n . 

We have not suceeded to get the answer to our question, but these results suggests 
that the question will be positive. 

4-1. Two-dimensional cases 

To construct low discrepancy sequences, we need a transformation not only expand- 
ing but also shuffling. For two-dimensional cases, let sq be the infinite sequence of 
0's, and 




Y,\C w \r^ <oo. 



W 




91 = 



WlW 2 • • • , 
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where w n is the word with length 2™~ 1 and only the last symbol equals 1 and other 
symbols are 0, that is, 

si = 101000100000001 

We consider the digitwise sum modulo 2 on the set of infinite sequences of and 1. 
Let 9 be a shift map to left. Then 

8 ,8i,9si, . . . ,9 n ~ 1 si 

generate all the words with length n. Let us define 

where we identify x G [0, 1) and its binary expansion x\Xi • ■ • . Then we can prove 
for I = [a, a + 2~ n+nl ] x [(3, [3 + 2~ n ~ m ) for binary rationals a and (3 and m < n,the 
image of I by F n does not overlap and F n (I) = [0, l) 2 . For example, I = [0, 1) x 
[0, 1/4) (n = m = 1), though 9y belongs only to [0, 1/2), x\ takes both and 1. 
Hence, 9y expands all [0, 1) by adding so or s± depending on where x belongs. Thus 
F(/) = [0,1) 2 . 

Theorem 4.1 ([7]). The van der Corput sequence generated by the above F is of 
low discrepancy. 



4-2. Three-dimensional cases 

For three-dimensional cases, we have not yet proved to get low discrepancy se- 
quences by one transformation. So we use a sequece of transformations F, F%, F%, 

Even for higher dimensional cases, we will be able to construct low discrepancy se- 
quences in the same way. However, we have no proof yet. 
We denote 

f - (!) - (? 

Then expressing x' , y' , z' in binary expansions, we define 



V'l 



■'-2 

V'l 

-,1 



V'J 



{Xn+l\ 

Vn+1 
Zn+1 



X n +2 
Vn+2 
Zn+2 



M 



Vi 
z, 



Xl 

yi 

Z\ 



Here, M is an infinite-dimensional matrix. One example of M is following. 
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X— 1 


V— 1 


Z — 1 


X — 2 


V — 2 


Z-2 


X — 3 


V — 3 


Z-3 


a:_4 


V — 4 


Z_4 


X\ 





1 


1 





1 








1 











1 


yi 


1 


{) 


1 








1 








1 


1 


{) 





Z\ 


1 


1 





1 








1 











1 





X2 





1 











1 





1 








1 





V2 








1 


1 














1 








1 


2*2 


1 











1 





1 








1 








x 3 





1 








1 











1 





1 





yz 








1 








1 


1 














1 


z 3 


1 








1 











1 





1 








x& 

















1 





1 











1 


y& 











1 














1 


1 








z& 














1 





1 











1 





X$ 














1 








1 








1 


1 


ys 

















1 








1 


1 





1 


Z5 











1 








1 








1 


1 





Xq 























1 











1 


VQ 


























1 


1 








Zq 




















1 











1 





x 7 























1 








1 





V7 
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*7 
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X8 
































1 





v& 
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28 





























1 








Xg 
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V9 
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Z9 
































1 





xio 






































yio 





{) 







































































This M has following properties: for any nonnegative integers n and m: the 
determinant of the minor matrix with coordinates 

{^l, ■ ■ ■ , ^n+m} X {j/— l, ■ ■ ■ , y—rn Z—l ■ ■ ■ i Z— mj 

and the determinant of the minor matrix with coordinates 

{xi , . . . , x n , yi , . . . , y m } x {z— ij . . . 5 Z— n — m } 

do not vanish. M also has symmetry in the permutations x, y and z. This matrix 
corresponds to sq and s± in two-dimensional case, and this shuffles coordinates, and 
we get for k > n + m 

I = [a, a + 2- fe -™-" 1 ) x [/?, P + 2- fe+n ) x [ 7 , 7 + 2- fe+m ) 

or 

I = [a, a + 2- k - n ) X [/3,/3 + 2- fc -' m ) x [ 7 , 7 + 2- k+n+m ), 

then 

^O0 = [o,i) 3 , 

where a, (3 and 7 are binary rationals. 

Theorem 4.2. The van der Corput sequence generated by F%,F2,... is of low 
discrepancy. 
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